The gravitational-wave memory from eccentric binaries 
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The nonlinear gravitational-wave memory causes a time-varying but nonoscillatory correction to 
the gravitational-wave polarizations. It arises from gravitational waves that are sourced by gravita- 
tional waves. Previous considerations of the nonlinear memory effect have focused on quasicircular 
binaries. Here, I consider the nonlinear memory from Newtonian orbits with arbitrary eccentricity. 
Expressions for the waveform polarizations and spin-weighted spherical-harmonic modes are derived 
for elliptic, hyperbolic, parabolic, and radial orbits. In the hyperbolic, parabolic, and radial cases 
the nonlinear memory provides a 2.5 post-Newtonian (PN) correction to the leading-order wave- 
forms. This is in contrast to the elliptical and quasicircular cases, where the nonlinear memory 
corrects the waveform at leading (OPN) order. This difference in PN order arises from the fact that 
the memory builds up over a short "scattering" time scale in the hyperbolic case, as opposed to 
a much longer radiation-reaction time scale in the elliptical case. The nonlinear memory correc- 
tions presented here complete our knowledge of the leading-order (Peters-Mathews) waveforms for 
elliptical orbits. These calculations are also relevant for binaries with quasicircular orbits in the 
present epoch which had, in the past, large eccentricities. Because the nonlinear memory depends 
sensitively on the past evolution of a binary, I discuss the effect of this early-time eccentricity on 
the value of the late-time memory in nearly circularized binaries. I also discuss the observability of 
large "memory jumps" in a binary's past that could arise from its formation in a capture process. 
Lastly, I provide estimates of the signal-to-noise ratio of the linear and nonlinear memories from 
hyperbolic and parabolic binaries. 

PACS numbers: 04.25. Nx, 04.25.-g, 04.30.-w, 04.30.Db 



I. INTRODUCTION 

Gravitational waves (GWs) are usually thought of as 
purely oscillatory phenomena. For example, the GWs 
from the coalescence of compact-object binaries tend to 
have a characteristic structure: as a GW passes through 
a detector, the frequency and amplitude increase, but at 
late times, the amplitude exponentially decays to zero 
from its peak value. However, the GWs from a variety of 
sources display nonoscillatory components as well. The 
simplest example is the GWs produced by the scattering 
of two unbound masses in a hyperbolic orbit [lj . Other 
examples include the GW signal from the asymmetric 
ejection of matter or neutrinos d, 0] in supernova 
explosions or gamma-ray-bursts @, HI ■ The GW signals 
from these sources all exhibit a property called GW mem- 
ory. This refers to a long time scale difference in the val- 
ues of the observed metric perturbation associated with 
the GW: 
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where /i+. x are the two GW polarizations. In a GW 
detector that is truly freely-falling (one that follows 
geodesies), a GW with memory can cause a permanent 
deformation in the detector (hence the term memory). 
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The nonoscillatory sources mentioned above are all 
examples of the linear memory 0-[ll[. In these cases, 
the nonoscillatory component to the GW arises from 
nonoscillatory motions of the source that are encoded in 
the matter stress-energy tensor T^ v . Because unbound 
gravitating systems have sources that undergo nonoscil- 
latory motions, the linear memory tends to occur in such 
systems. In addition to the linear memory, there is also a 
nonlinear memory effect [T^l - fToT l . This nonlinear memory 
arises from the gravitational waves produced by gravita- 
tional waves: it is sourced not by nonoscillatory motions 
encoded in T^ u but rather in Tjjff [H,[I3; which describes 
the stress-energy of radiated GWs. These radiated GWs 
are themselves always "unbound" and hence produce a 
memory [l8( | . Furthermore, since this effect originates di- 
rectly from the radiated GWs and not the motion of the 
source, the nonlinear memory is present in all sources of 
GWs, including bounded systems. 

Why study the nonlinear memory? One of the key 
goals of gravitational-wave astronomy is to experimen- 
tally probe our understanding of general relativity (GR). 
While solar system and binary pulsar tests can only 
probe GR in the weak-field regime where the theory is 
nearly Newtonian, observations of coalescing compact- 
object binaries (and especially merging binary black 
holes) will produce GWs whose properties will depend 
on the strong-field, highly-dynamical sector of GR. There 
are a variety of nonlinear effects which will imprint them- 
selves on the waveforms from such systems. Of these 
effects the nonlinear memory is unique because (i) it 
describes waves produced by waves (among the most 
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nonlinear of interactions present in the theory), (ii) its 
nonoscillatory imprint on the GW signal is distinct from 
the manifestation of other nonlinear effects, and (hi) al- 
though it arises from a higher-order nonlinear interaction, 
it can enter the post-Newtonian (PN) expansion of the 
waveform amplitude at leading (OPN) order. 

The nonlinear memory also has the property of being 
hereditary. This means that the nonlinear memory piece 
of the waveform amplitude at some observer's time t de- 
pends not only on the configuration of the source at the 
corresponding retarded time t — R (where R is the dis- 
tance to the source) , but rather on the entire past-history 
of the source. To illustrate this property, compare the 
leading-order "quadrupole" expression for the GW held, 

hj? = ^lJ k T (T R ), (1.2) 
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with the corresponding expression for the nonlinear mem- 
oryd, 
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where ljk is the mass quadrupole moment, 
GW energy flux, rij is a unit radial vector, 



(1.3) 

-mm is the 
TV is a unit 

vector pointing from the source to the observer, and TT 
means to take the transverse-traceless projection (see also 
Sec. I of [13 )■ In Eq. (jl.2[) . one can clearly see that the 
primary GWs measured by a detector are directly deter- 
mined by the retarded time configuration of the source 
(as encoded by its mass quadrupole). But in the expres- 
sion for the nonlinear memory [Eq. (|1.3[) ]. its value at 
time Tr depends on a integral into the infinite pastQ So 
to determine the value of the nonlinear memory at one 
instant, one needs to know the energy flux (and hence 
the motion of the source) at all previous times. 



A. Motivation 

Previous calculations of the nonlinear memory have fo- 
cused almost exclusively on quasicircular binaries. Early 
work focused on computing the memory during the inspi- 
ral phase for quasicircular orbits, starting with Wiseman 
and Will [190 (see also 123]). Remarkably, they found 
that the nonlinear memory affects the GW polarizations 
at leading (OPN) order: 
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1 The GW tail effect exhibits a similar property, but in that case, 
the integrand drops off more steeply in the past than does the 
memory, so only the nearby past need be considered (see Sec. 4 
of [2l| or Sec. 5.3.4 of [22l | for a more detailed discussion). 

2 Note that the expression for the memory in the text after Eq. (17) 
of [Tgjl is missing a factor of 2, but the curves in their Fig. 1 agree 
with the expressions below. 
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where M = mi + mi is the binary's total mass, r\ — 
711x1712 /M 2 is its reduced mass ratio, x = (Afu;) 2 / 3 is 
the standard PN expansion variable for quasicircular bi- 
naries, (p(t) is the orbital phase, uj = <p is the orbital 
angular frequency, c@ = cos 8, sq = sinG, and (6,$) 
are the angles in the source frame that point in the ob- 
server's direction N. The second term in Eq. (|1.4a|) is the 
nonoscillatory memory term. (For a convenient and stan- 
dard choice of the polarization triad, there is no nonlin- 
ear memory contribution to the x polarization.) In [20(, 
these nonlinear memory contributions were computed to 
3PN ordeiQ, thus completing the PN expansion of the 
waveform amplitude consistently to that order. 

More recent work has investigated in detail the non- 
linear memory produced by merging binary black holes. 
Using a simple analytic model as well as an effective- 
one-body [24 1 approach, the full evolution of the memory 
for the inspiral, merger, and ringdown was computed in 
[25l [26| and the prospects for its detection with inter- 
ferometers were examined^ Calculations of the nonlin- 
ear memory from numerical relativity were performed in 
[27L [28j , and its detectability with pulsar-timing-arrays 
was considered in [29j - [3l| . 

The purpose of this study is to analyze the behavior 
of the nonlinear memory for arbitrarily eccentric bina- 
ries. There are several reasons why this generalization 
is worth considering. First (and most importantly), real 
binaries will be eccentric and not quasicircular. Even 
though gravitational radiation-reaction tends to circular- 
ize binaries, we expect to observe some binaries with non- 
negligible ellipticity (see, e.g., Section I and Appendix A 
of [32[ as well as [3J, |34j)' or with hyperbolic trajecto- 
ries (from scattering events in clusters or galactic nuclei) . 
However, even if one is only interested in nearly circular 
binaries, the hereditary nature of the nonlinear mem- 
ory makes it important to consider the eccentric case as 
well. Because binaries which are currently quasicircular 
were more eccentric in the past, this prior eccentricity 
has modified the orbital motion and could potentially in- 
fluence the calculation of the nonlinear memory integral 
in Eq. (| 1 . 3[) . Clearly, if we hope to actually observe the 
nonlinear memory effect, we must also be prepared to 
account for binary eccentricity. 

A second motivation comes simply from the desire to 
have a complete and consistent understanding of the 
waveforms produced in the general two-body problem. 
For example, waveform polarizations for elliptical bina- 
ries are implicit in the classic work by Peters and Math- 
ews [HI (although they focus on computing the radi- 
ated power), and are given explicitly first in Wahlquist 



3 Arun et. al [2~D | previously showed that the 0.5PN memory con- 
tribution vanishes. 

4 Thorne also examined the memory's detectability, treating 
it as an unmodeled burst, while Kennefick |2.'j| considered only 
the inspiral contribution to the memory. 
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[36j and later in Refs. |37H39| to leading (OPN) order. 
These elliptical waveforms have since been extended to 
1PN order in amplitude by Junker and Schafer [IS] (see 
also 0, EH), to 1.5PN order in Blanchet and Schafer 
[HI, and to 2PN order (neglecting tails) in Gopakumar 
an d Iyer Frequency-domain waveforms are given 

in [32|, I45l448| . (The GW phasing for elliptical binaries 
is currentl y k nown to relative 3PN order in the conser- 
vative [49HiI| and dissipative parts 0, 0, 51, H H- 
6l|.) As we will see later, even the leading-order 
(Peters-Mathews/Wahlquist) waveforms are incomplete 
because — as in the quasicircular case — the nonlinear 
memory modifies the polarizations at OPN order. 

In the case of hyperbolic orbits, waveform polariza- 
tions were first derived at OPN order in Turner [l|, with 
1PN corrections computed in (42[. In the large- 
eccentricity (bremsstrahlung) limit, waveforms were com- 
puted to 1PN order in [55|,[62| and using the "post- linear" 
formalism in (63l . [64j . These waveforms already show a 
linear memory, but the nonlinear memory has only been 
calculated in the large-eccentricity limit [19| . Waveforms 
for radial orbits are also considered up to 1PN order in 
551 (although the energy flux has been computed at 2PN 
6a j and 3PN (66j j orders), but the nonlinear memory con- 
tribution in the radial case has not yet been computed. 

In this paper, I will attempt to complete our knowledge 
of the nonlinear memory for binaries with arbitrary ec- 
centricity, including elliptical, hyperbolic, parabolic and 
radial orbits. Unlike in Ref. [2(| where the nonlinear 
memory was computed to 3PN order in the quasicircular 
case, in this paper I will restrict the calculation to the 
leading-PN-order piece of the nonlinear memory. I will 
also discuss how the nonlinear memory (because of its 
hereditary nature) is affected by the prior eccentricity of 
a nearly circularized binary. 



B. Summary 



In Sec. Ill A[ we begin by reviewing the prescription in 
.20f] for computing the nonlinear memory from the GW 
energy flux. This involves decomposing the GW polar- 
izations into a sum of spin-weighted spherical harmonic 
modes hi m7 and relating the nonlinear memory modes 
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(mom) 
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to "lower-order" oscillatory modes hf that are 



accurate to Newtonian (OPN) order. These Newtonian- 
order modes depend only on the mass quadrupole mo- 
ment 1-xm , for which explicit expressions are easily derived 
for general Newtonian binaries (Sec. Ill Bj) . In Sec. Ill BT1 
formulas for Keplerian orbits are reviewed, and general 



5 Additional works also consider the GW polarizations in the case 
of eccentric binaries with spinning components, but for simplic- 
ity, I will not consider spin effects here. I will also not discuss 
results based on black hole perturbation theory, except to say 
that the nonlinear memory has not been considered in that for- 
malism. 



expressions for Keplerian waveforms are derived. The 
material in this section is concisely presented and may 
be useful to those interested in leading-order waveforms 
valid for any eccentricity. 

In Sections III CI III Dl and III El I specialize the nonlin- 
ear memory calculation to the cases of elliptical, hyper- 
bolic, parabolic, and radial orbits. The primary results 
are explicit expressions for the nonlinear memory modes 
and the corresponding polarizations /i+, x - Aside 



(mom) 



from these explicit expressions, the following are some of 
the primary results of this analysis: In the case of inspi- 
ralling elliptical binaries, the nonlinear memory behaves 
similarly to the quasicircular case. The primary contri- 
butions come from the m = modes, which have the 
scaling 
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where pit) is the semilatus rectum of the ellipse and Ti m 
is a hypergeometric function that depends on the eccen- 
tricity et(t) and its value e_ at some early time [see Ap- 
pendix [B] for the exact expressions, or Eqs. (12.35)) for the 
low-eccentricity limit]. Note that as in the quasicircu- 
lar case, the nonlinear memory modifies the waveform 
at leading (Newtonian) order [cf. Equation (ll.4a|) with 
M/p ~ 0{x)\. As discussed below, this result would 
hold even if we were to extend our analysis to orbits that 
undergo periastron advance. In the case of hyperbolic 
and parabolic orbits, this scaling is very different: 
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Hi m (e t ), (1.6) 



where Him [ which can be read off of Eqs. (|2.51|) ] is a 
function of the eccentricity et (which is constant in the 
absence of radiation reaction) . Note that in this case 
all of the (l,m) modes contribute to the nonlinear mem- 
ory, which is a factor of rj(RI/p) 5 / 2 smaller than in the 
elliptical case. (A similar scaling also holds in the case 
of radial orbits.) 

It is easy to understand the reason for this difference 
in scalings. The nonlinear memory can be written as a 
time integral of the form 
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where the time derivative h 



dhi^^/dt has 



the same leading-order scaling [cx ry 2 (M/p) 5 ] for all or- 
bits. The difference in the scalings results from the time 
scale on which the integration is carried out. In the case 
of elliptical orbits, one must integrate over the entire in- 
spiral, which occurs on a radiation-reaction time scale 
T„, so that 
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In the hyperbolic case, one effectively integrates over a 
much shorter "scattering" or "orbital" time scale T or b0 



A/i (mem),hyperb. _ ^-XD^ _ h fa«*X» M (j£ 
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Since 
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(1.9) 
(1.10) 



we can see why the nonlinear memory in the hyper- 
bolic/parabolic case enters at a much higher PN order 
than the elliptical/quasicircular case. Integrating over 
the much longer radiation-reaction time effectively allows 
the memory to "build up" to a much larger value than 
one would naively expect from such a high-order PN ef- 
fect. 

One of the more important results from this study con- 
cerns the analysis of the dependence of the memory on 
the early-time history of the binary (Sec. IIIip . For ex- 
ample, all quasicircular binaries presumably had ellipti- 
cal orbits earlier in their evolution, and this increasing 
(into the past) ellipticity could eventually result in a hy- 
perbolic binary. To address this issue, I computed the 
evolution of the dominant mode of the nonlinear mem- 
ory waveform as a function of time for a nearly circu- 
lar binary. This evolution was computed (i) assuming 
that the binary always remains quasicircular, and (ii) as- 
suming that the binary's eccentricity increases as time 
evolves into the past. A comparison of these two evo- 
lutions is shown in Fig. [S] (where time is parameterized 
in terms of the eccentricity of the elliptical binary). Ac- 
counting for the evolving eccentricity makes a small (but 
non-negligible) correction to the memory. The eccentric- 
ity correction is small because even though eccentricity 
is increasing into the past, so is the orbital separation (or 
the semilatus rectum). Since the integrand in Eq. (II. 7[) 
is weighted by a factor of (M/p) 5 , its value drops off at 
larger separations; so late-time values (smaller p, when 
et is also small) are weighted more heavily than the dis- 
tant past (when the eccentricity is large, but M/p is very 
small) . 

An eccentric binary could also produce a large memory 
jump in its distant past, for example due to its sudden 
formation in a capture process. Such a memory jump 
would in principle be observable, but only if one's GW 
detector is operating when (in retarded time) that jump 
occurred. In general, effects of the early-time history of 



To be precise, one is still integrating over the infinite time that 
the hyperbolic orbit spans, but the contribution to the memory 
integral mostly builds up over the short amount of time it takes 
the binary to "scatter" or significantly change its direction. This 
"scattering time" is not precisely defined, but it scales the same 
way (via Kepler's law) as what we would normally call an "orbital 
time" in the case of bound orbits. 



the binary are not observable in the memory signal if they 
occurred before the start of the observation. Hereditary 
effects are therefore only important over the observation 
time, and one need not worry about knowing the state of 
the binary prior to the start of the observation. This is 
illustrated with an explicit example in Sec. lIIII and Fig.|6] 
Lastly, in Sec. IIV1 I discuss how to use the waveforms 
for hyperbolic obits to make rough estimates for the 
signal-to-noise ratio (SNR) of a memory signal from a 
gravitational-scattering event. Unlike the case of bound, 
inspiralling binaries (where the innermost stable circular 
orbit or the ringdown provides a natural high-frequency 
cutoff), the SNRs from hyperbolic binaries depend on a 
distance of closest approach (as well as an eccentricity 
parameter) and are dominated by the linear rather than 
the nonlinear memory. For future ground-based detec- 
tors, linear and nonlinear memory signals from the scat- 
tering of stellar-mass binaries within our galaxy are po- 
tentially detectable. Space-based detectors could see lin- 
ear memory signals from the scattering of supermassive 
black holes (SMBHs), but detecting the nonlinear mem- 
ory from these events will be more difficult. (This is in 
contrast to the nonlinear memory from merging SMBHs, 
which should be detectable to moderate redshifts 

Some useful results are presented in the Appendices. 
Appendix lAl derives general expressions for the spherical 
harmonic modes of the mass and current source multipolc 
moments at Newtonian order. Appendix [B] shows how 
to express certain integrals over eccentricity in terms of 
hypergeometric functions. Appendix [C] discusses the role 
of wavelength averaging of the GW stress-energy tensor 
in the nonlinear memory calculation. 



II. COMPUTING THE NONLINEAR MEMORY 

We begin this section by first reviewing the decom- 
position of the GW polarizations in terms of multipole 
modes hi m on a spin-weighted spherical harmonic basis 
fSec. HI A| . The nonlinear memory pieces of these modes 



(mem)-] 



dtdil, 



are related to integrals over the GW energy flux 
This is described in [2(| , and the reader is directed 
there for a more detailed exposition. Since we are only 
interested in the leading-order memory, it is sufficient to 
express the energy flux in terms of the Newtonian-order 
nonmemor 30 modes hf m . Section |HB] gives explicit ex- 
pressions for these Newtonian-order modes and special- 
izes them to Keplerian orbits parameterized in terms of 
the true anomaly angle v. These modes are then sub- 
stituted into the expression for the nonlinear memory 
modes (Sec. iHBll) . Sections |HCl |HDl and HH] special- 
ize the calculation to the case of elliptical, hyperbolic, 



7 When referring to "nonmemory" modes here and below, I really 
mean the "non-nonlinear-memory" pieces of the modes, which 
could contain oscillatory or linear memory pieces. 



5 



parabolic, and radial orbits, providing explicit expres- 
sions for the /i[™° m ^ 
and x polarizations. 



sions for the hiT, e modes and their corresponding 



A. General expressions for the waveform and 
nonlinear-memory modes 

The GW polarizations can be decomposed as a sum 
over multipole modes via 



h + -ih x =2_^ 



/i lm _ 2 Y lm (e,$) , 



1=2 m=- 



where 



G 



V2Rc l + 2 



U tm (T R ) - -V m (T R ) 
c 



(2.1) 



(2.2) 



Here, _2F' m (0,$) are the spin- weighted spherical har- 
monics, (R, O, $) are the distance and angles that point 
from the source to the observer, and U lm and V lm are the 
spherical harmonic representations of the radiative mass 
and current multipole moments (see Sec. II A of [20] for 
details and notation). Constructing the GW polariza- 
tions requires explicit expressions for the -2Y lm (Q, $). 
The general formula for the _ 2 i^ m (0, $) can be found 
in Eqs. (2.13) and (2.14) of [20|. Here, we will only need 
the following modes: 



20 



>Y 



2±2 



'-(l±c e ) 2 e 



iY*" = - 



y4±2 _ - 




±2i* 



1), 



y4±4 _ 



'-(l±c e )'(lT7c e + 7c4)e 



-s^(l±ce) 2 e ±4l$ , 



±2i* 



(2.3a) 
(2.3b) 
(2.3c) 
(2.3d) 
(2.3e) 



where we define ce = cosG and Se = sinG. 

At leading order, the radiative moments Ui m and Vi m 
reduce to the source moments J; m and Ji m . At higher PN 
orders, the radiative moments are corrected by tail terms 
and other nonlinear couplings [see, e.g., Eqs. (2.21), 
(2.22), and (2.32) of (2(|]. Since we are focused on only 
the leading-order memory contribution, we will ignore 
all of these higher-order terms except for the nonlinear 
memory contribution C//™ em ' ) itself. Note that there is no 
nonlinear memory contribution to Vi m . For our purposes, 
we can therefore ignore all current multipole moments. 
This allows us to approximate the waveform modes at0 

htoK^c+hSr). (2.4) 



Here, the nonlinear memory piece is given by [see 
Eqs. (2.32) and (3.3) of 20]] 



*-'-TrVp+^y^*y dn 3 (n)1B,<n) 



16tt /(/-2)! [ Tr 



dE„ 



m-\-m 



A Kj Vl"lm'-m"-r 



R dt(h llm ,h* l/)m)) ), (2.5) 



where in the second line we have substituted the expan- 
sion for the e nerg y flux -j^Pj in terms of the hi m modes 
[Eq. (2.28) of M] and G^ 2 % 3 mi m 2 m 3 k an angular inte- 
gral proportional to the product of three spin- weighted 
spherical harmonics (see Appendix A of |20|). The an- 
gle brackets ( ) mean to average over several wavelengths 
of the GW and arise from the averaging implicit in the 
construction of a well-defined GW stress-energy tensor 
[l6l [I?} (see Appendix [C] for a discussion of the implica- 
tions of this averaging) . 

Note that the nonlinear memory modes /i[™ om ^ are 
themselves defined in terms of the full hi m modes. But 
in practice, the "nonlinear memory contribution to the 
nonlinear memory" is negligible, so we only need sub- 
stitute the nonmemory modes into the right-hand-side 
of Eq. (|2.5p . Furthermore, since we are only concerned 
with the leading-order nonlinear memory, we need only 
substitute the leading-PN-order piece of the I = 2 mode 
/i2m into the right-hand-side of Eq. (|2 . 5[) . 

Now we proceed to compute the angular integrals in 
Eq. (|2.5|) and explicitly express the in terms of 

the individual /i 2m modes. We begin by noting that the 
angular integration implies certain selection rules on the 
maximum I for which the are nonzero (see Sec. Ill 

B of [ITJ]). Since our calculation is to leading-order, only 
the /i2m modes enter the right-hand-side of Eq. (|2.5p . The 
selection rules then imply that only /4m 

, (mem) 



4r m> > and 



will be nonzero; all higher- 1 nonlinear memory 
modes vanish (this was checked by explicit calculation). 

Our results are further simplified by the fact that the 
nonmemory piece of the hi m modes are approximated by 



'Hm 



r(l) 



(2.> 



where the N emphasizes that our results are valid only 
at Newtonian order. Since If = (— l) m /;_ m , the non- 
memory modes also satisfy hy n * = (— l) m /iJL~,- We also 



note that since 



Ira V / I — i n' 

.- K ^/m(^' f) ( see Appendix [A]), and if 
we assume that the orbit lies in the x-y plane (so that 
9 = ir/2), then /i^ti ^ ^2±i ^ sin 20 = 0. These simpli- 
fications imply 



We denote n th time derivatives by A( n ' = d n A/dt r ' 



;N .N* 
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<H^2 1 "2±1 



. N* 
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0. (2.7) 
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Defining h 



mcm)(l) 



dhi™ em ' I (ITr, explicitly evaluat- 
ing the angular integrals in Eq. (|2.5p (with hi m — > hf m 
on the right-hand-side) then yields 
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(mem)(l) 
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(mem)(l) 
4 40 



(mcm)(l) 
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(mem)(l) 
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= h 



(mem) (1) 
3m 



= h 



(mem) (1) 
4±1 



= h 
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1260 V 2tt 
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2o2 V 2 ^^ o/l 2W 



i? /14 
504 V 2^ 



(//. 



N \2 
2±2J 



= o, 

(2.8a) 
(2.8b) 

(2.8c) 

(2.8d) 

(2.8e) 

(2.8f) 



m 2 = *J 2 ivMr 



2-2 M r, 

r ip =p zirrip 

r I 



r(3) _ 



r( 3 ) - 4 

J 2±2 — ~ 4 
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2tt /M 



5 ^ \ r 



(f =p Atrip) . 



)• 

(2. lid) 
(2. lie) 

(2.11f) 



These expressions are valid for general orbits that satisfy 
the Newtonian equations of motion (|2.10[) . The explicit 
Newtonian-order polarizations are found from substitut- 
ing Eqs. (|2.6p and (|2.3|) into Eq. (|2.ip and summing only 
over the I = 2 modes n 



(i + 4) 



•2 2-2 

r — r tp 



M 



cos2(ip — $) 



-2rr<p sin 2(ip - $)] - Sq ( r 2 + r 2 ip 2 



t) 



(2.12a) 



B. Explicit expressions for the hf m and frj™™ 
modes for Newtonian binaries 

Now we write out explicit expressions for the hf m 
modes. A derivation of the source mass and current mul- 
tipole moments for Newtonian binaries is given in Ap- 
pendix [A] The result for the mass quadrupolc moment is 
found in Eq. (IA15al) . The I = 2 moments for an orbit in 
the x-y plane are 



(2.9a) 



4 N ±2 = 2yy V Mr(t) 2 e^\ (2.9b) 

where M = mi + m2, r\ = mxmijM 2 , r(i) is the relative 
orbital separation, and ip(t) is the relative orbital phase. 

Next we compute time derivatives of these mass mo- 
ments. We eliminate second derivatives using the New- 
tonian equations of motion: 



2 M 
r = rp j, 

r 2 



2rip 



(2.10a) 
(2.10b) 



Dropping the "N" label, the resulting derivatives of the 
moments are 



h« = 2^c e 



R 



2frip cos 2(ip — $) 



•2 2-2 

r — r ip 



sin 2(1/? - $) 



(2.12b) 



Substituting h$ m = I^/{R^/2) and Eqs. (j2"TTj) into 
Eqs. (|2T5)1 then gives 



• (mem)(l) 
'20 



, (mem)(l) 
"2±2 



21i?V 15 



16 ^(£)V + 24rV)), 



(2.13a) 



21RV 5 ' \ V r 



(2.13b) 



h 



(mem) (1) 
40 



(2.13c) 



, (mcm)(l) 
"4±2 



(2.13d) 



(mem) (1) 
4±4 



551/7"' ((t > ^ 4 "<™ 2 



(2.13e) 



r(i) 



7T 

— r/Mrr, 
15 



(2.11a) 

7 2±2 = A \J^n Mr (* T ir<f) e T2iv , (2.11b) 



(2.11c) 



for the nonvanishing memory modes. 



9 These formulas agree with Eqs. (6) of 1531 if we choose <& = ±7r/2 
and change the overall sign on both polarizations. This difference 
arises from the choice of the polarization tensors (see Sec. IIA of 
[20] for the convention used here). 



FIG. 1. (color online). Notation and parameters describing elliptical and hyperbolic orbits. The left-hand figure shows a 
particle moving on an ellipse at a distance r from the focus and making an angle v (the true anomaly) with respect to the 
pericenter. The ellipse has an eccentricity et and a size described by either its semimajor axis a, focus-pericenter distance r p , 
or semilatus rectum p. Also indicated is the eccentric anomaly u, which is the angle from the pericenter to the projection of the 
particle's position on the circle that circumscribes the ellipse. The right-hand figure shows a particle moving on a hyperbolic 
orbit, where we additionally indicate the asymptotes of the hyperbola (dashed lines), the particle's impact parameter b, and 
the scattering angle Q s . The argument of pericenter w is taken to be zero in these diagrams. 



1. Formulas for Keplerian orbits 

Now we wish to specialize these expressions to Kep- 
lerian orbits. The time evolution of the orbital separa- 
tion r and phase tp is parameterized in terms of the true 
anomaly v (not to be confused with the orbital speed V) , 



P 



(2.14a) 
(2.14b) 
(2.14c) 



where, for planar orbits, only three orbital elements are 
needed to parameterize the binary: the semilatus rectum 
p, the eccentricity^ et, and the argument of pericenter 




10 Throughout this article, we denote the eccentricity by et to avoid 
confusion with the mathematical constant e and to emphasize 
that the eccentricity can evolve with time pes., et = et(t)]. This 
choice is not meant to imply an identification of our eccentric- 
ity parameter with the "time eccentricity" used in the quasi- 
Keplerian formalism that describes PN elliptical orbits (e.g., [52j 
and references therein). In that formalism, three eccentricity 
parameters, et, e r , and e v , are introduced; but at Newtonian or- 
der these three eccentricities are equivalent and we can identify 
either of them with the et used here. 



zd. Figure Q] illustrates the meaning of the various orbital 
parameters introduced throughout this article. Time is 
determined by integrating Eq. (|2.14c[) : 



t-U 



dv 



M J (1 + e t cosv) 2 
The following derivatives follow from Eqs. (|2.14[) : 



(2.15) 




ip-r), 



2M 



-e t suit). 



(2.16a) 

(2.16b) 
(2.16c) 



Note also that the instantaneous orbital velocity is given 
by 



2,u2 



V =r + r z ip- 



M 
P 



(1 + e 2 t + 2e t cosv), (2.17) 



the orbital energy per reduced mass is 
E V 2 M IM 



2 p 



(1-e?), 



(2.18) 
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and the semilatus rectum is related to the semi-major 
axis a by 



p = a\l - e\\ 
and to the pericenter distance r p by 
P = r p {l + e t ). 



(2.19) 



(2.20) 



Applying Eqs. ff2TT4]t and ([2T6]l . Eqs. (j2~TT|) simplify 



to 



r(2) 
z 20 



7T M 2 . 

—X] e t (e t + cosv), 

15 p 



r(2) 
l 2±2 



(2.21a) 



x [1 - e 2 + (1 + e t cosu)(l + 2646*™)] , (2.21b) 



r(3) 
£ 2±2 



2tt /M\ 



5/2 



p) 



(1 + e t cost;)^ 



X [e t sin v =F + e t cos «)] . (2.22b) 

The modes that appear in the waveform polarizations 
[Eqs. (|2.2ip ] are plotted in Fig. [5] for hyperbolic and 
elliptic orbits. To produce these plots, the differential 
equation for the true anomaly [ (12. 14c[) ] was solved nu- 
mericallvF^I Note that the hyperbolic waveforms show a 
linear memory in the (2, ±2) mode. This is discussed in 
more detail in Sec. Ill PI 

To compute the polarizations in terms of the true 
anomaly, one substitutes the following expressions into 
Eqs. (12321 : 

M 

frip= — (1 + et cosv)et sinw, (2.23a) 
P 

MM, . , 

— = — e t (e t +cosv), (2.23b) 
r p 

- = (2 + 3e t cosD + e 2 cos2w). 

r p 

(2.23c) 



rV 



r ip 



r(3) 
'20 



it /M\ 5/2 . N2 

— rj [ — e t smu 1 + e t cosv) 
15 V P J 



(2.22a) 



Likewise, we can use the above expressions for r, r, and 
Lp to write Eqs. (|2 . 1 3|) in terms of e*, p, and v: 



(mem)(l) 



32 

21R V 30 



— J (1 + e t cos w) 4 (24 + e 2 + 48e t cosu + 23e 2 cos 2 u) ), 



l 2±2 



(mcm)(l) 
40 



16 
2LR 

2 /2tt 



2LRV 5^ 



315i? 



^1 P T^ 
P 

M 



et sin u(l + et cos i>) 4 [et sin v =F 4i(l + et cos i>)] 



ft 



(mcm)(l) 
4±2 



315i? 



M 



(l + e t cosv) (16 + 3e t + 32e t cost; + 13e t cos v) ), 



— I e T2lv e t sinu(l + e t cosw) 4 [e t sinw =p 4i(l + e t cosw)] ), 



(2.24a) 
(2.24b) 
(2.24c) 
(2.24d) 



'4±4 



-^Wy)] 2 / f — j e T4iv (l + e t cosv) 4 [16 - e t 2 + 32e t cosv + 17e t 2 cos 2 w ± 8ie t sinw(l + e t cosu) 



(2.24e) 



where the angle brackets again arise from the averaging 
inherent in the definition of the GW energy flux. 



11 These waveforms could also be produced by using a represen- 
tation in terms of the eccentric anomaly u. In this case, the 
problem of solving for the time involves finding a root of the so- 
called Kepler equation (rather than solving a differential equa- 
tion). However, in contrast to the true anomaly parameteriza- 
tion, separate sets of equations must be used to treat the elliptic, 
hyperbolic, and parabolic cases (see, e.g., Ch. 6 of [H). 



C. Elliptical orbits 

For bound, eccentric orbits (0 < e t < 1), the wave- 
length averaging in Eqs. (|2.24j) is accomplished by ex- 
plicitly averaging over an orbital period P rb- For any 
function F(t), this orbit-averaging is defined by 



(F(t)) 



1 



dtF(t) 



orb Jo 

(1 - e 2 ) 3 / 2 > 2T 
2tt 



(2.25) 



dv 



F(v) 



(1 + e t cost;) 2 ' 
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FIG. 2. (color online). Waveform modes for hyperbolic and eccentric orbits as a function of a normalized time coordinate 
i — (t/M)(M/p) 3 ^ 2 . Both plots show the Ii m modes [Eqs. (|2,21a[) and (|2.21b|) ] that appear in the construction of the leading- 
order waveform polarizations. The eccentricity is indicated in each plot and completely determines the shape of the curves. I 
have assumed that the binary passes through the pericenter angle zn = at t = 0. Note that the imaginary part of the I22 mode 
shows a linear memory in the hyperbolic case (the size of this memory scales like et/p ~ 1/6 for large et and impact parameter 
6; see Sec. Ill D")) . Parabolic orbits (e t = 0) produce waveforms very similar to the hyperbolic ones, but all modes asymptote to 
the i-axis at t — > ±00 (no linear memory). The 72-2 modes are similar to those plotted and are given by Re72-2(£) = Re/22 (i) 
and Im/2-2(i) = — Im/22(— t). 



where 



Pn 



2tt 



(2.26) 



(1 - e 2 ) 3 / 2 V M 
follows from Eq. (|2.15[) . Averaging Eqs. (|2.24l) then yields 



. (mcm)(l) 
l 20 



(mcm)(l)_ 
L 2±2 



(mem)(l)_ 
40 



256 HT 2 /Mf f1 2 , 3 / 2 /; 145 2 73 I 

(2.27a) 

(2.27b) 

^0 (1 -e?,-( 1 + le ?+ ^ei) 

(2.27c) 

,(mc m )(l)_ 13y^ 2 /M 5 2 2 3/2 / 2 2 \ 2m 

(2.27d) 

%5 



52 k 
21R 

64 
315i? ' 



(mom)(l)_ 
M±4 — 



^A/7^ l ~ )ef(1 ~ e * )3/2eT4iro ' (2 ' 27e) 



This averaging has essentially removed the high fre- 
quency (<~ 1/Porb) structure from the memory waveform. 
Appendix [C] explicitly shows that this averaging proce- 
dure has only a very small effect on the memory. 

Next, we need to compute the time integrals of the 
above expressions. In the circular case treated in pfj| |. 
this was accomplished by changing variables to x = 



(Mw) 2 / 3 . In the eccentric case, both the eccentricity et 
and the semilatus rectum p vary with time, but p can 
be easily expressed in terms of et- We therefore change 
variables from time t to eccentricity, and integrate from 
some early-time value of the eccentricity e_ to its value 
at some later time e + = e t (t): 



• (mem) 
l lm 



, (mem)(l) 
'Hm ai 



■+ h 



(mcm)(l) 
Im 

det/dt 



-det. (2.28) 



In computing the above integral, we will need to make 
use of the following equations for the evolution of p and 
et [these are easily derived from the results of [H, [6!| , 
along with the Newtonian-order relations in Eqs. (|2.26[) 
and ([2~T9]) ]: 

| = -^(y) 3 (l-e?) 3 / 2 (8 + 7 e ?), (2.29a) 



rj ( M 



de, 



dt 15M V P 



e t (l - e t 2 ) 3/2 (304 + 121e 2 t ). 



Dividing the first equation by the second yields 
dp _ U P (8 + 7e t 2 ) 



(2.29b) 



(2.30) 



de t 'e t (304+ 121e 2 )' 
which can be solved to give 

p( et ) = E± et WW (304 + me 2 ) 870 / 2299 , (2.31) 
Co 
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where 



C = e 12/19 (304+ 121eg 



2\870/2299 



(2.32) 



and po is the value of p at some arbitrary reference time 
when et = eo. Since both p and et evolve with time, this 
relation allows us to eliminate the time-dependent p(t) 
terms in Eq. ()2.28|) and instead express the integrand 
entirely in terms of the evolving eccentricity e* (which is 
our new integration variable) and the constants po and 
eo- The integrand in Eq. (|2.28p for the relevant modes 
then becomes: 



. (mom) 
l 40 



5n V M 2 

56 Rp * 



1 



7T r\M 2 



63 V 10 Rpo 



ir rjM 2 



1 



63 V 10 Rp 



1- — 



12/19 



26/19' 



1-1^- 



£o 

e_ 

12/19' 



(2.35b) 



12/19' 



(2.35c) 



dh 



(mem) 



det 



2 /IOtt rjM 2 C (192 + 580e 2 + 73e 4 ) 



7 V 3 R Po e 3Vi9 (304 + 121e 2 ) 3169 / 



2299 ' 

(2.33a) 



dh 



(mom) 
2±2 



de, 



dh 



(mem) 
40 



4^5^ rfM 2 C e[ /19 (13 + 2e 2 ) 

7~ Rp (304 + 121e 2 ) 3169 / 2299 ' 

(2.33b) 

1 fTr/M 2 Co (128 + 396e 2 + 51e 4 ) 



rfe, 



<2/i 



(mem) 
4±2 



(mem) 
4±4 



de t 



42 V 10 i?po e f /19 (304+ 121e 2 ) 3169 / 2299 ' 

(2.33c) 

(2.33d) 

(2.33e) 



V?7?M 2 C e[ /19 (13 + 2e 2 ) 
21 Rpo (304+ 121e 2 ) 3169 / 2299 

25 fnriM 2 



O e t 



24 V 7 R Po (304+ 121e 2 ) 3169 / 2299 



These expressions can be analytically integrated and the 
result expressed in terms of hypergeometric functions. 
Details of this are given in Appendix |Bl where we show 
that all of the above modes can be expressed in the form 



F Im (e_)], (2.34) 



= A lm C (e )e* im ™ [F lm (e t ) 



where Ai m are constants that can be read off of Eqs. (|B1[) , 
and Fi m is a sum of hypergeometric functions given in 
Eq. (|B3[) . Note that in computing the integral over e* 
we have chosen the integration constant such that the 
memory vanishes at an early-time eccentricity value of 
e_. We have also ignored periastron precession (choosing 
n7 to be fixed), but the relaxation of this assumption will 
be discussed below. 

In the limit of small et, we can easily evaluate the 
integrals of Eqs. 



. (mem) 
l 20 



2 / IOtt tjM 2 
7V 3 Rpo 



12/19 



12/19' 



2 jlO^TjM 2 

7V ~3 Rp 



(mem) 
l 2±2 



v/Stt f]M 2 
56 Rpo 



e e 



e_ 
eo 



12/19' 



26/19 



(2.35a) 



26/19" 



(mem) 
l 4±2 



672 i?p 



e^e 2 



26/19 



672 i?p 



1-1 — 

e* 



e_ 
eo 

26/19" 



26/19" 



l 4±4 



25^ rjM 2 T4lro 4 
172 032 i?p 6 ° 

25^77? V M 2 4 - ro 4 
172 032 i?p 6 6 * 



64/19 



1 - I — 

et 



e_ 
eo 

64/19" 



(2.35d) 



64/19" 



(2.35e) 



e 12/19 304870/2 299 j and in thc 



where we have used Co 
second line of each equation we have re-expressed the re- 
sult in terms of the time-dependent p(t) « Po(et/eo) 12 / 19 
using Eq. (|2"3Tj) . 



In the zero-eccentricity limit, ^±2" 



■ (mem) 
1 4±2 



and 



(mem) 



^ 4 ™ 4 m ' ) vanish and h^ em ^ and /i 4 ™ em ' ) reduce to the 
circular-orbit values found in Eqs. (4.1)-(4.3) of [2(j. In 
the elliptic case, the m ^ modes do not contribute to 
the memory for several reasons: first, they are suppressed 
by factors of e™ or eg", and their numerical coefficients 



-0 1 

tend to be much smaller than the coefficient of h 



(mem) 
20 



More importantly, the factor of e T "" ro is not actually 
constant as we have assumed so far. Post-Newtonian cor- 
rections result in periastron precession, which causes w 
to vary with time [or with changing et(t)\. For example, 
the rate of periastron advance is 



2?rfc 



orb 



i(f) S V e ^ (2.36) 



for k = 3M/p at leading-PN-order. Using Eqs. (I2.29b|) 
and (|2.31[) . the pericenter angle can be obtained as a 
function of eccentricity by solving 



dm 45 / po \ 
de t ~ r) VM/ 



3/2 



°o e t 



(304 + 121e 2 ) 994 / 2299 ■ 



(2.37) 



For arbitrary (bound) eccentricities, this equation can 
be integrated using Eq. (|B2I) . For small eccentricity, this 
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FIG. 3. (color online). The left-hand plot shows the ft^o" and memory modes as a function of eccentricity e t for two 

values of the reference eccentricity eo and for e_ = 1. The dashed lines show the small-eccentricity approximation [Eqs. (|2,35[) ]. 
In practice, these curves would not extend to arbitrarily small et but would terminate when the last-stable-orbit is reached. 
The right-hand plot shows the evolution of the h ™ mode for various values of the indicated reference eccentricity eo (every 
curve passes through the point et = eo when r p = 20Af). The integration is terminated at the last-stable-orbit (at which point 
the curves flatten). In a real merger, the memory would continue to evolve past this point, eventually saturating at a different 
value. 



equation has the solution 



W = ZUq 



_5_ /po\3/2 

32?7 \ MJ 



18/19' 



(2.38) 



where vuq is the value of w when e t = eo- As the eccen- 
tricity varies from e to 0, the pericenter angle changes 
by Am ~ (l/rf)(po/M) 3 ' 2 ; this leads to oscillations in 
Eqs. (|2.33[) that, upon integration, cause the m ^ terms 
to be further suppressed. Hence, as in the quasicircular 
case, only the m = terms contribute to a secularly in- 
creasing memory effect. 



In the left-hand plot of Fig. |3l we plot the h 



and h 



(mem) 
40 



(mem) 
20 



modes as a function of e*. This is obtained 



from both the full analytic solution for the mode evo- 
lution [Eq. (|2.34l) . solid lines] and the low-eccentricity 
limit [Eqs. (|2.35l) . dashed lines], choosing e_ = 1 in both 

mode is much smaller than 



cases. Note that the h^ em 
the h^"^ mode. 



It is also convenient to express the above results 
in terms of the pericenter distance r p rather than p. 
The time evolution of r p is found from differentiating 
Eq. (j2~20| and using Eqs. (j2~29| : 



dr-r. 



7] /M\ 3 (l-e t ) 3/2 



(192-112e t +168e 2 +47e 3 ). 



dt 15\r p J (1 + et) 7 / 2 ' 

(2.39) 

The evolution with eccentricity r p (et) is easily found from 
Eq. (H3U), 



>'o 



12/19 



C (l + e t ) 



(304 + 121e 2 ) 870/2299 (2.40a) 



12/19 
e 

(l + e ) 



(304+ 121e 2 ) 870/2299 . 



(2.40b) 



This allows Eqs. ([2~2"T]) . (|2~g3l . and ([235)) to be expressed 
in terms of r p or tq. 

The right-hand plot of Fig. [3] attempts to further il- 
lustrate the dependence of the h^ " 1 ^ memory mode for 
different eccentricities. In place of a time variable t, we 
can parameterize the temporal evolution in terms of the 
eccentricity et. This is because, at Newtonian order in 
the conservative dynamics, an inspir ailing eccentric bi- 
nary passes at some point in its evolution through every 
value of et S (0, 1) with a one-to-one mapping between 
t and et (assuming we neglect the details of the binary's 
formation or its interactions with the external universe). 
To distinguish one eccentric binary from another, we need 
to specify the value of the eccentricity at some fiducial 
orbital separation. The different curves in the right-hand 
plot of Fig. [3] are parameterized by the value of the eccen- 
tricity e t = eo when the binary passes through a pericen- 
ter distance of r p = 20 M. The curves are obtained from 
the analytic solution for the (2,0) mode in Eq. (|2.34l) . 
choosing e_ = 1. The (2,0) mode is allowed to grow 
until the last-stable-orbit (LSO) is reached, correspond- 
ing to the condition p = r p (l + e t ) = 6 + 2et [7fJ [the 
LSO value of eccentricity eLSO is determined by com- 
bining this condition with Eq. (|2.40j) ]. This plot shows 
that binaries with large eo reach the LSO while they are 
still mildly eccentric and at slightly smaller values of the 
memory. (However, note that GWs radiated during the 
merger and ringdown will cause the memory to grow past 
its LSO value.) 
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The polarizations for the nonlinear-memory waves for 
bound, eccentric orbits can be simply computed by sum- 
ming the m = modes in Eq. (|2.1j) : 



(mem) 



(mem) 



1 /30 

8 



(mem) 



V_3 .(mem),, 2 

n 40 [icq 



= 0. 



(2.41a) 
(2.41b) 



For arbitrary eccentricities, Eq. (|2.34|) for the I = 2 and 
I — 4 modes must be substituted into the above equation. 
In the small-eccentricity case, Eqs. (|2.35|) yield 



^+ ncm) = ^4(17 + 4) 



48i?p 
IV I(17 + C|) 



12/19 



e_ 

12/19 



12/19 



(2.42) 



Note that in the circular limit this agrees with 
Eqs. (4.4)-(4.6a) of J2C| or Eq. ([i~4a)l above. 



D. Hyperbolic and parabolic orbits 

To treat the case of hyperbolic and parabolic orbits, 
we ignore the possibility of periastron advance and we 
fix the periastron direction to lie along the +x axis. In 
this case, the reduced mass particle swings around the 
origin in a counterclockwise sense, entering at very early 
times along the asymptote at ip = v_ = — arccos(— l/e t ), 
and exiting at very late times along the asymptote at 
if = v + = arccos(— l/e t ) (see Fig. [1]). The corresponding 
scattering angle S is given by 



2arccos(— l/e t ) 



(2.43) 



For hyperbolic orbits, it is also useful to define two ad- 
ditional parameters that can be used in place of p or e t ■ 
The asymptotic velocity is 



V 2 



E M 



P 



(el -I). 



(2.44) 



The impact parameter 6, defined to be the perpendicular 
distance from the center of mass M to the ingoing or 
outgoing asymptote of the hyperbola, is found to be 



6 = 



P 



1 



(2.45) 



The equations below can alternatively be expressed in 
terms of Voo or b using the above relations. 

It is important to note that the waveforms from hy- 
perbolic orbits already contain a linear memory [l|. For 



example, consider Eqs. (|2.21al) and (|2.21b[) for et > 1 and 
ip = v varying between V- and w+ . This difference in the 

orbital phase angle does not affect the 1^ mode (since 
cosu is even), but it does affect the imaginary part of 

(2) 

^2±2 (which is odd), leading to a memory in that mode 
(see Fig. . More explicitly, the linear memory jump be- 
tween late and early times for a hyperbolic orbit is found 

(2) (2) 

from the difference I 2 m{v+) ~ -^2m( w -)i yielding [see also 
Eqs. (10) of H3] 



Ah 



(lin. mem) 



Ah 



(lin. mem) 
2±2 



o, 



±il6w — r) 



7T MM(ef-l) 3 / 2 



(2.46a) 
(2.46b) 



5 ' R p ef 
The corresponding memory jump in the polarizations is 



p, mH= MM(^ (i + 4) 



R p 



Ah ,nn. mem; = _ 4 „-- V"« ^ (l + ci)sin2$, 

(2.47a) 

a , (lin. mem) , M M (eg - if' 2 

Ah\ = —on— 5 cecos2$. (2.47b) 

R p e$ 

Note that for large et0 



M (ef - if 2 
P e 2 t 



M 



(2.4 



and Eqs. (j2T47|) agree with Eqs. (15) of jl9|. Note also 
that for parabolic orbits (et = 1), the linear memory 
vanishes. This is because the asymptotic incoming and 
outgoing directions of the orbit are now the same (u_ = 

— 7T, V + = +7r). 

To compute the nonlinear memory we proceed from 
Eqs. (|2.24p . Since unbound orbits are no longer peri- 
odic, there is no need to average over an orbital pe- 
riod. Instead, we directly perform the time integrals over 
Eq. (|2.24[) . changing variables to the true anomaly using 
Eq. (|2"3gc|) : 



. (mem) 
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v(Tr) t(mem)(l) 

dv. (2.49) 



The integrand is a sum over powers of sines and cosines 
or their products, and is easily evaluated for any limit 
of integration. For simplicity (and because periastron 
passage happens relatively quickly for hyperbolic orbits), 
we will focus on computing only the overall memory jump 
Ah[™ cm \ rather than the evolution of the memory with 



time: 
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The resulting memory modes for any et > 1 are 
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Note that there is a memory contribution of order r\ 2 {M jp) 1 1 2 in each mode (a relative 2.5PN correction). Note also 
that each mode is real- valued. These modes are plotted in Fig. |4j The resulting polarizations are 
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For parabolic orbits (et — 1), Eqs. (|2.5ip simplify to: 
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(2.53a) 

(2.53b) 
(2.53c) 
(2.53d) 
(2.53e) 



12 In this case, we have ~ e 2 (M/p) and 6 ~ p/et; or alterna- 
tively, e< « V£(6/M) and p « V^^/M). The e t > 1 limit 
then corresponds to the bremsstrahlung (small-angle scattering) 
limit, Vjjj, 3> M/b. Note that the scattering angle for et S> 1 is 
e s » 2/e< n (2/V£)(M/6) « 1. 



and the corresponding polarizations are 
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(2.54b) 

Unlike in the linear-memory case, the nonlinear memory 
for parabolic orbits is nonzero. Even though parabolic 
orbits are marginally bound, their radiated GWs are un- 
bound and hence contribute to the nonlinear memory. 

We can also examine the e t 3> 1 limit. In this case, it 
is easy to extract the large-et behavior from Eqs. (|2.51l) : 
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FIG. 4. (color online). The A/i;™ em ' modes for hyperbolic or- 
bits from Eqs. (|2.51|l . For e t = 1, these reduce to Eqs. (12.531) . 
while for e t 3> 1, they asymptote to Eqs. (|2.55|) . 



a > nonlin.mem 2 ^ I ^ I y 



,M 

~R 



(2.59) 



This indicates that the nonlinear memory for high- 
velocity gravitational scattering is typically much smaller 
than the linear memory (see also Sec. IIV|) . This is in con- 
trast to the case of bound eccentric (and circular) orbits, 
where the linear memory vanishes (but see Sec. V B of 
20] for a caveat) while the nonlinear memory is oc r\. 
These scaling differences in the nonlinear memory arise 
from differences in the integration time over which the 
nonlinear memory builds up [see the discussion following 
Eq. (fTTT)) above]. 



E. Radial orbits 

Next we consider radial orbits corresponding to the 
head-on collision or separation of two masses. In this 
case, the equations of motion and conserved energy yield 



ip = (p = 0, 
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(2.60a) 
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the polarizations are given by 
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(2.55a) 
(2.55b) 
(2.55c) 
(2.55d) 
(2.55e) 

(2.56) 
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[192 cos 2$ + s|(1756 - 128 cos 2$ - 50 cos 4$) 

-s|(102 - 32cos2$-25cos4$)] , (2.57a) 



A , (mem) 7T 2 M ( M\ 

x [96 sin 2$ - s| (16 sin 2$ + 25 sin 4$)] . (2.57b) 

This agrees exactly with Eq. (16) of [19| . providing fur- 
ther confirmation of the correctness of the above results. 

Note also the different scalings between the linear 
[Eq. (|2.46[) ] and nonlinear memories in the et S> 1 limit: 

M M 

(2.58) 
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(2.60b) 



The multipole modes in Eqs. (|2.1ip easily simplify in the 
radial case (where we can choose ip — const = 0), and 
the leading-order waveform polarizations become 

hN + = v -f{( r2 -v) [d + 4) cos 2$ -4]}, 

(2.61a) 
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If the relative radial velocity approaches Woo at infinite 



separation, r — M/ 



r -> v" 



+ M/r. Radial waveforms 



can therefore show a linear memory effect that depends 
on Voo and the initial and final values of M/r. 

To compute the nonlinear memory, we simplify 
Eqs. (|2.13p (again choosing ip = 0). We easily see that 
all of the leading-order memory modes have the form 
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Clm R It 
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where the constants C/ m can be read off of Eqs. (12.131) . 

Converting the time-integral of /il™ " 1 ^ 1 ' to a radial in- 
tegral and using Eq. (|2.60bl) . 
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= ±V2jE + M/r, 



(2.63) 



(where the upper sign here and below indicates increasing 
radial separation and lower sign indicates radial infall), 

the h& modes can be expressed as 



M 
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E+—dr, (2.64) 
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where r± refers to the value of r at late or early times. 
Evaluating the integral yields 
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For the case of radial infall from rest at infinity, the 
above simplifies to 
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where r + — > r(t). We can also consider the case of ra- 
dial binary disruption (e.g., this could also model a star 
that radially ejects a piece of material). If the initial 
separation is r_ and the asymptotic late-time relative 
velocity of the two components approaches Woo (so that 
E = u^/2), the resulting nonlinear memory shift is 
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Note that in both of the above cases the nonlinear mem- 
ory is a relative 2.5PN correction to the Newtonian wave- 
form. In all cases, the waveform polarizations for radial 
orbits are given explicitly by 
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" :cm) is given by Eqs. (I2J51) . (I2J61) . or (12^671) 
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mode. (For simplicity and because they tend 
to be much smaller, I will neglect the other memory 

modes.) In the first model, we consider the ft,2o 1Gm '° lllp '' 
mode for an elliptical binary described via Eq. (|2.34[) . 
with e_ = 0.99 and eo = 0.01 at a pericenter distance of 
ro = po/(l + eo) = 6M. This mode is plotted as the blue 
(solid) line in Fig. [SJ We will also need to model how 
the eccentricity evolves with time. To do this, I evolve 
Eqs. (|2.29b|) and (|2.39[) . but I change to a new time vari- 
able T = — t so that I can more easily evolve the system 
"backwards" in time starting from the initial conditions 
r p (T = 0) = r = 6M and e t (T = 0) = e = 0.01. For an 
equal-mass binary, I find that the "early-time" eccentric- 
ity e_ = 0.99 is reached at a time T_/M w 2.031 x 10 8 . 
This mode is compared with a purely quasicircular model 



for the h. 
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mode which is given by 
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where 
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(3.2) 

This model forces both 



and r_ = r(T_) 225.8M. 
the quasicircular mode /i2Q lom ' clrc ^ and the elliptical mode 
/i 2 ™ em ' olllp '' to vanish at the same time (T = XL), which 
can be considered the start of the observation. It also 
ensures that both orbits have a pericenter separation 
r p = ro = 6M at time T = 0. The two modes are 
plotted in Fig. [SJ where the value of time for both modes 
is parameterized in terms of the eccentricity of the ellip- 
tical mode. This figure indicates that the quasicircular 
model provides a moderately accurate representation of 
the true evolution of the memory mode (which accounts 
for the orbit's past eccentricity). At the end of the evolu- 
tion (T = 0, r p = 6M), the two modes have a fractional 
error of s» 1.5%. 

To better quantify the degree to which the two modes 
"overlap," I have computed the following normalized in- 
ner product: 



O = 



J n T_ hi{t)h2(t)db 



Chj(t)dt] \£-hj(t)dt 



(3.3) 



III. SENSITIVITY OF THE MEMORY TO THE 
EARLY-TIME HISTORY OF A BINARY 

In this section, we wish to evaluate the degree to 
which the nonlinear memory from a quasicircular inspi- 
ralling binary is sensitive to its deviations from circu- 
larity. These deviations arise from the binary's initial 
eccentricity, and are damped by radiation reaction (in 
absence of external perturbing forces). To perform this 
evaluation, we compare two models for the evolution of 



where h\ = h 



(mcm,circ) 



and ho = h 



(mcm.cllip) 



This is 



equivalent to the commonly computed overlap between 
two GW signals, but here assuming white noise. For val- 
ues of e = 0.01 or 0.001, I find the value O « 0.976; this 
decreases slightly to 0.975 for eo = 0.1. Although I have 
not considered a realistic noise model, this calculation 
suggests that ignoring the effects of past eccentricity in 
quasicircular binaries is a reasonable approximation and 
is not likely to result in significant reduction in the signal- 
to-noise ratio of the nonlinear memory. 
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FIG. 5. (color online). Evolution of the (2,0) memory mode 
in the quasicircular (red/dashed) and elliptical (blue/solid) 
cases for an equal-mass binary. The two modes are evolved 
as described in the main text. Time for both curves is param- 
eterized by the eccentricity of the elliptic case. The elliptic 
mode is set to zero at an eccentricity of e_ = 0.99 and evolves 
until it reaches a value of eo = 0.01 at r p = 6M. The quasi- 
circular mode also evolves up to a separation of 6M and is set 
to zero at the time when the eccentricity of the elliptic mode 
is 0.99. Two of the insets zoom in on the early- and late- 
time stages of the evolution (at high and low eccentricities, 
respectively) . The other inset shows the absolute value of the 



difference between the two curves, \h. 



(mem.circ) 



■ (mcm,cllip) 



Now let us consider the memory that results over the 
entire lifetime of a binary system, including its initial for- 
mation. As one would intuitively expect, any bound ellip- 
tic binary experiencing gravitational radiation-reaction 
evolves to larger eccentricities (and larger orbital separa- 
tions) into the past until e< > 1 and the binary becomes 
unbound. This was proved rigorously in (tTJ . Equiva- 
lently, for certain choices of its initial orbital parameters, 
a hyperbolic binary can lose energy from gravitational- 
wave emission and become bound. The waveform for 
such a scenario can be approximately modeled using 
Eqs. (|2.21al) and (|2.21bl) combined with a prescription 
for the instantaneous evolution of the orbital elements 
(sec, e.g., [71]). If we choose our time and angular co- 
ordinates conveniently so that capture happens at peri- 
astron, a schematic description of the waveform modes 
from such a captured binary would look like the left- 
hand plot of Fig. d] for i = (t/M){M/pf/ 2 < smoothly 
matched onto the right-hand plot of Fig. [2] for t > 0. 
(Note that the different modes and their slopes in that 
figure have the correct signs at t = to allow for such a 
matching.) After capture, such a binary would circular- 
ize and eventually merge, with the waveforms evolving in 
the standard way for t > 0. Of course, this description 
is somewhat idealized. In the real world, other inter- 
actions (e.g., tidal dissipation, three-body interactions, 
gas drag, or dynamical friction) are more likely to result 



in binary capture (although gravitational radiation losses 
could play an important role in very dense stellar systems 
such as globular clusters or galactic nuclei) However, 
for the purpose of considering the size of the memory 
jump over very long time scales, let us presume that at 
some early time the binary is in an unbound, hyperbolic 
orbit, while at some later time it is a bound, elliptic bi- 
nary that circularizes and merges. For such a binary, the 
total memory jump is roughly given by Eq. (|l.ll) with 

riM 2 (e 2 - l) 3 / 2 

^^ k VV"' (3 - 4a) 

AE 

lim h+.x oc — — , (3.4b) 

t->oo Jri 

where pi and are the semilatus rectum and eccentricity 
prior to captur J^l. and AE is the energy radiated in GWs 
throughout the inspiral, merger, and ringdownj^fl This 
suggests that large memory jumps can result not only 
from the nonlinear memory (which grows most rapidly 
during the final phases of coalescence), but also through 
the linear memory associated with binary capture. 

A more relevant issue is the observability of some sig- 
nature of the formation or early-time state of the binary. 
Clearly, if one's GW detector is operating when the bi- 
nary capture process occurs (in retarded time), then the 
signature of the capture, including the resulting memory, 
will be seen by the detector (provided it is sensitive to 
low- frequency effects like the memory). However, what 
if the capture process (and the associated passing GWs) 
occurred long before the start of the observation period? 
Does the capture process still leave an "imprint" on the 
waves observed at later times? Intuitively, one expects 
the answer to this question to be "no." This is indeed 
correct as can be seen with the following argument. 

Consider a simplified GW detector consisting of two 
particles floating in space separated by a distance L. 
Placing the first particle at the origin of its own proper 
reference frame, the position x 3 of the second particle 



13 Another possibility is that the binary was not captured but was 
"born bound," with each component star forming from a frag- 
menting molecular cloud. The system could then have evolved 
into a compact-object binary and a source of GWs. 

14 The eccentricity after capture is approximately given by eo = 
e* + Ae, where Ae = -i^r?(M/pi) 5 / 2 . This can be derived 
by considering the change in eccentricity for a parabolic binary 
(e< = 1), 

(3.5) 



Ae : 



/7T 
-7T 



dv 



-dv, 



where det/dv = (det / dt) / (dv / dt) . An expression for the instan- 
taneous (i.e., not orbit-averaged) value of det/dt can be derived 
by considering the Lagrange planetary equation (cf. 68]) for an 
osculating Keplerian ellipse under the action of the 2.5PN radia- 
tion reaction force (see also the last equation in [7l| or Eq. (2.14) 
of[U). 

However, note that the nonlinear memory is only proportional to 
the radiated energy at leading-order in an (2, m) mode expansion 
of the energy flux |2(J . 
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FIG. 6. (color online). Response of a freely falling 
gravitational-wave detector to a memory signal. The differ- 
ential displacement [normalized by the detector length L and 
asymptotic value of the memory h}- + °°^] is plotted vs time. 
The memory sig nal is modeled by Eq. ((37TT|) with /i (_oo) = 
and t = 10M. The second time-derivative in Eq. (|3.12[) is also 
plotted (blue/dotted curve). The solid (black) curve shows 
the response of a detector that observes the entire memory 
signal [cf. Equation (|3.9[l ]. The dashed (red) curve shows the 
response of a detector that starts monitoring the memory at 
time to — — lOr (indicated by the arrow) [Eq. (|3.10l) ]. In this 
case the t < — lOr build-up of the memory is lost, and a linear 
drift (proportional to the slope at to = — 10r) develops at late 
times. 

relative to the first is given by the equation of motion 
(Ch. 35.5 of 0) 



-h 



TT fc 

jk x ) 



(3.6) 



where overdots here refer to the derivative of the proper 
time at the first particle, and hJjF is the metric perturba- 
tion in transverse-traceless gauge. We can choose to ori- 
ent our two particles and the resulting coordinate system 
such that their motion along their direction of separation 
x is given by 



x{t) = ^h+(t)x(t). 



(3.7) 



For very small displacements, x(t) « L + Sx(t), and the 
equation for the difference in the particles' relative sepa- 
ration simplifies to 



5x(t) = -h+(t). 



(3.8) 



Now we consider two scenarios: in the first we assume 
that our detector has been freely-floating for all times, 
so it observes the entire build up of the memory. In 
the distant past, we assume that our memory signal ap- 
proaches the value and that its derivative vanishes, 



-oo) — > 0. In this case, Eq. (|3.8I) has the solution 
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Already we see in this case that the value of the GW field 
in the asymptotic past [h^ 00 ^] is not observable; instead 
only the difference between that asymptotic value and 
the current value (at time t) is observable. However, since 
the detector has been operating for arbitrarily long times, 
the measured value of the memory retains any imprint of 
the past evolution of the binary [e.g., the value of the 
nonlinear memory at time t would depend on the motion 
of the source all the way to t — > — oo, but not on the 
value of 

In the second scenario, let us suppose that the memory 
signal has started arriving, but our detector is rigidly 
fixed in position until some time to when we allow our 
particles to_be free-floating. In this case, Eq. (13.81) has 
the solutior 
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Sx(t) = | 



(mem) . 



(3.10) 

Here (somewhat obviously) we see that the memory loses 
its dependence on times before t . We also see the de- 
velopment of a linear drift (proportional to the slope of 
the waveform at to)- This drift arises from the initial im- 
pulse the detector receives from the passing wave at the 
moment it is released. 

To make the above discussion more explicit, consider a 
schematic model for a nonlinear memory waveform given 
by the arctangent function: 



(mem) 



h (+co) _ h (-oo) / t \ h (+co) + h (-oo) 

arctan^-j + , 

where h^°°' are the asymptotic values of the memory 
and r is the characteristic rise time of the memory. The 
second time-derivative of this function is 



■(mem) = 2[/ I (+°°> - hS^jtr 
l + TT^+T 2 ) 2 



(3.12) 



These functions and the resulting differential displace- 
ments are plotted in Fig. [6] for the two scenarios men- 
tioned above. This graphically illustrates that we can 
only observe the build up of the memory that occurs 
while our detector is operating. A similar model could 
also be based on the hyperbolic tangent function, 
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2 ' 
(3.13) 



16 This situation is equivalent to solving Eq. 1 13. 8 6 with the right- 
hand side multiplied by a Heaviside function Oft — to). 
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which approaches its asymptotes more quickly and has a 
second derivative given by 



;(H = [fa (+oo) - fe^ 00 '] smhft/r) 
l+ t 2 cosh 3 (i/r) 



(3.14) 



IV. 



ESTIMATING THE SIGNAL-TO-NOISE 
RATIO OF MEMORY JUMPS 



Here I provide some simple formulas for estimating the 
detectability of the memory. For the case of merging qua- 
sicircular binaries, detectability estimates are presented 
in 126] and will be discussed in more detail in 27]. For 
elliptical binaries, we have seen that the memory behaves 
quite similarly to the quasicircular case, so the estimates 
of detectability are little changed. I instead will focus 
on the linear and nonlinear memory for hyperbolic and 
parabolic orbits. 

We begin by defining the angle-averaged square of the 
signal-to-noise ratio as 



(P 2 ) = 



h 2 c (f) df 
hlif) f ' 



(4.1) 



TABLE I. Evaluation of the quantity N in Eq. (|4~6|) for dif- 
ferent detectors. The integral is taken from a low-frequency 
cutoff of 30 Hz in case of initial LIGO, 10 Hz for advanced 
LIGO and advanced Virgo, 1 Hz for the three iterations of the 
Einstein Telescope (ET), and 3 x 10~ 6 Hz for LISA. In the 
A?max column, the upper limit of the integral [f c in Eq. (|4.6[l ] 
is taken as 2000 Hz for the ground-based detectors and 1 Hz 
for LISA. In the iVi ow column, the upper limit of the integral 
is taken to be the frequency where fh% is minimized. The 
sensitivity curves used here can be found in [78 
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transform is given b 



\h+Af)\ 



Ah, 



2tt/ 



(4.3) 



However, a real memory signal has some finite rise time 
r which imposes a hi gh- frequency cutoff at f c ~ 1/r in 
the Fourier transform£j We can therefore approximate 
the characteristic strain by 



-(|AM 2 + |AM 2 } 1/2 e(iW). 



(4.4) 



where the average is over all sky positions, source 
orientations, and polarization angles [see, e.g., 
Eqs. (2.33)-(2.36) of [z3]. Here, h n (f) = y/afS n (f) is 
the sky-averaged rms noise amplitude per logarithmic 
frequency interval. The factor a is 5 for orthogonal 
arm detectors like LIGO and 20/3 = 5/ sin 2 (60°) for 
equilateral triangles like LISA or the Einstein Telescope 
(73} . The characteristic amplitude is given bsFl 



/a/) = 2/(iM/)r + iM/)r) 



2\l/2 



(4.2) 



where a tilde denotes a Fourier transform. If we approx- 
imate the memory as a step-function, then its Fourier 



17 For cosmological sources, one must replace / — > (1 + z)f and 
R — » Dl,(z)/(1 + z) in this expression, where z is the redshift 
and Dl is the luminosity distance. 



This follows from the Fourier transform of the Heaviside function, 



f 

J —I 



' H(±t)e 2m f t dt= S -^-±-L. 
K ' 2 2tt/ 



This step-function approximation is equivalent to the zero- 
frequency limit (ZFL) discussed in Rj |74h76|| . In that case, one 
approximates the Fourier transform of the time-derivative of a 
signal h(t) near / ~ via 

/•-(-co r-\-oo 

h(f) = / h(t)e 27rift dt / h(t)dt, 

J OO J OO 

= h(+oo) - h(-oo) = Ah, 

and we use_ the usual relation for the Fourier transform of a 
derivative, h(f) = (— 2nif)h(f), to arrive at 

iAh 



2tt/ 



19 For an explicit example of this, consider the Fourier transform 
of the signal in Eq. ll3J3t for f > [77t: 

r(mem) [h<+°°> - h<-°°>] . . . a n 

h + = «7TTCSCh(7T TJ) 

iffc(+oo) _ J,(-°a)l f _4 -1 

= - J f ' {i-^^ + oK^]}. 

Here, one can see from the Taylor expansion the sharp cutoff in 
the ZFL value of the Fourier transform when / ~ 1/r. 
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The SNR then becomes 



(4.5) 



where we define h c to be Eq. (|4.4|) without the Heaviside 
factor and 



N 



df_ 

fK 



-1/2 



(4.6) 



To evaluate N, one needs to choose a value for the cut- 
off frequency f c ~ r _1 . The rise time for a hyperbolic 
trajectory is r = nM(p/M) 3 ' 2 , where K is a factor that 
depends on how the rise time is defined 23 Although N 
clearly depends on the parameters of the binary (through 
f c ), we tabulate its value for several detectors in Table U 
by choosing f c to be either the location of the minimum 
of fh 2 n or the high-frequency cutoff for the detector. 

For the case of a pulsar timing array (PTA), we can 
estimate N by making use of Eq. (31) in [30j : 



A, 



PTA 



„ lfi /250 20 
1.94 x 1(T 16 — — 



1/2 /lOyrs 



(- 



V 100ns/ 

where A t is the number of measured timing residuals, N a 
is the number of pulsars in the array, a n is the noise in 
the timing residuals (assumed to be Gaussian stationary 
white noise that is uncorrelated and the same for each 
pulsar), and T Q b s is the total observation time for the 
PTA. The numbers used above are for near-future PTAs. 
The Square Kilometre Array (SKA) [83j will achieve bet- 
ter sensitivity, including a factor of 10 decrease in the 
noise, and perhaps a factor of ~ 100 increase in the num- 
ber of suitable pulsars [8J|. 

Taking the angle-averages of Eqs. (|2.47l) . (|2.54[) . and 
(|2.57|) . and expressing the results in terms of the peri- 
center distance, yields the following characteristic ampli- 
tudes for parabolic and hyperbolic orbits: 



h 
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(4.8a) 
(4.8b) 
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v'21 (IT- - ) 91,. , ,, , ., , , 

h c,e t ^i = 5^ V — [ — ) v oo, (4.8c) 
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where the second and third equations show the lin- 
ear and nonlinear memory for hyperbolic orbits in the 
large-eccentricity limit (for which the memory is largest). 



20 If we define the rise time by taking the integral in Eq. 112.151 1 over 
v £ [— G s /2, s /2], then we find that n < 4/3 and asymptotically 



approaches 2/ef for et 2> 1. 



Note that in the hyperbolic case, the nonlinear mem- 
ory is smaller than the linear memory by a factor w 
0.79r)(M/r p ) 2 Voo- In practice, this amounts to a factor 
3> 4 orders of magnitude, so we ignore the nonlinear 
memory in the hyperbolic case. Plugging in numbers for 
some plausible (but perhaps optimistic) sources yields 



v (mem) 
l c,e t =0 



2.2 x 10~ 22 



r] \ 2 / M /IOMq \ /20M 
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(4.9a) 



— -"(i^(StX^). 

= 9.6 X 10-'» m g«\ 
V0.25/ V fl/lGpc J\ r p J 

(4.9b) 

[For cosmological distances we should take M/R — > 
M z /Dl in the above expressions, where M z = (1 + z)M 
is the redshifted mass.] The SNR can then be esti- 
mated by combining Eq. (|4.9j) with Eq. (|4.5|) and the 
numbers in Table HI These rough estimates indicate that 
GW bursts with linear memory could be detectable 
with second-generation ground-based detectors and fu- 
ture space-based detectors. The nonlinear memory from 
GW bursts from unbound (or marginally) bound bina- 
ries will be more difficult to detect and will likely require 
third-generation detectors. Current and near-term PTAs 
are not sufficiently sensitive to detect memory bursts of 
the types considered here. 



V. CONCLUSIONS 

This work has generalized previous computations of 
the nonlinear gravitational-wave memory effect to the 
case of binaries with arbitrary eccentricity. In the case 
of hyperbolic, parabolic, and radial orbits, the nonlinear 
memory is a 2.5PN correction to the waveform. In the 
case of elliptical binaries, the nonlinear memory affects 
the waveform at leading order (just as in the quasicir- 
cular case). To completely describe elliptical waveforms 
at leading-PN-order, the nonlinear memory contributions 
derived here should be added to the well-known nonmem- 
ory expressions first derived by Peters and Mathews [35[ 
and Wahlquist 36]. I have also investigated the sensi- 
tivity of the nonlinear memory to the early-time history 
of the binary. In the case of quasicircular binaries that 
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were initially elliptical, the early-time eccentricity pro- 
vides only a small correction to the memory. Further- 
more, I have shown that contributions to the memory 
made outside of the observation time are undetectable. 
Lastly, I provided simple estimates of the signal-to-noise 
ratio for memory bursts arising from sources on unbound 
orbits. 

There are a variety of areas in which this study could 
be extended. The nonlinear-memory calculations pre- 
sented here are restricted to leading order. For hyper- 
bolic, parabolic, and radial orbits the waveforms are only 
known to 1PN order, so there is little motivation to com- 
pute higher-order corrections to the leading-order non- 
linear memory terms (which themselves enter as 2.5PN- 
order corrections to the waveform). However, in the el- 
liptical case, the oscillatory waveform polarizations are 
known to 2PN order, so 2PN-order corrections to the 
leading-order nonlinear memory terms would be needed 
to have complete 2PN-order elliptic waveforms. In ad- 
dition, the effects of spinning binary components on the 
nonlinear memory have not yet been computed. This cal- 
culation is in progress in the case of quasicircular binaries 
and will be reported elsewhere [85j . Computing the non- 
linear memory for eccentric, spinning binaries will be left 
for future work. It would also be interesting to investi- 
gate the size of the linear and nonlinear memory in the 
case of ultrarelativistic collisions and scatterings [86U9l| . 
These situations could show a very large memory effect. 
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Appendix A: DERIVATION OF THE 
LEADING-ORDER MASS AND CURRENT 
SOURCE MOMENTS FOR A GENERAL 
TWO-BODY SYSTEM 

The purpose of this appendix is to derive expressions 
for the mass and current multipole moments in the form 
of (l,m) modes that are valid for any two-body orbit 
at Newtonian order. At leading order, we are only con- 
cerned about the so-called source moments which are de- 
fined in terms of integrals over a stress-energy pseudoten- 
sor. General expressions (valid for any PN order) for 
the mass and current symmetric-trace-free (STF) source 
multipoles, II and Jl, can be found in Eq. (85) of f92j . 



These STF tensors with L indices (where L = a\a^ ■ ■ ■ a{) 
can be difficult to work with, and for some calculations 
it is more convenient to instead use the "scalarized" ver- 
sions of these moments, Jj m and J; m . These "scalar" 
multipoles are simply the coefficients of the expansion 
of the STF mass and current multipoles on the basis of 
the STF spherical harmonics 3^ m [these are defined in 
Eq. (2.12) of [93[ and are related to the standard scalar 
spherical harmonics via Eq. (|A8[> below]. The STF mo- 
ments and their (I, m) modes are related by the following 
formulas [Eqs. (4.6) and (4.7) of [H]: 




and h m 

Jim 



Bij L y l i 



lm * where 



At = 



B, 



16tt + 2) 



(2/ + l)!!y 21(1-1) ' 
32ttZ / 



(I + 2) 



(2Z + 1)!! V 2/(Z + !)(/-!) 



(A2a) 
(A2b) 



(A3a) 
(A3b) 



Now we specialize the general form for STF mass and 
current moments in Eq. (85) of 92] to the Newtonian- 
order moments for general orbits (but arbitrary Z-value). 
This derivation could be easily extended to the 1PN- 
order moments (see [94| HI])- The Newtonian-order 
source mass and current multipole moments for a sys- 
tem of N (nonspinning) point masses is 



N 

E 

A=l 
N 



m A y A 



(L) 



Jl — mA£ Va v a 



(A4a) 



(A4b) 



A=l 



where A labels the body, the multi-index L refers to a 
product of I vectors (e.g., y\ — y^yT ' ' ' vV)i £abc ls ^ ne 
Levi-Civita tensor, and the angled brackets () mean to 
take the STF projection on the enclosed indices. The 
"N" superscript emphasizes that these are Newtonian- 
order moments. We now specialize to a two-body system 
with masses m\ and m.2, total mass M = m\ 4- and 
reduced mass ratio t] = toito 2 /M 2 . We transform to the 
center-of-mass frame using 



mi _ 

Vl = ~M X> 

mi _ 



V2 



M ' 



(A5a) 
(A5b) 
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where x = y\ — yi — rn has length r and n points from 
77i 2 to tyi\. We also define the individual and relative 
velocity vectors via va = Va and v — x. 

Substituting the above relations into Eqs. (|A4[) gives 
[Eqs. (5.21) and (5.22) of [H]: 



j 

rN 



l - rjMsi(r))x {l 



Jl = V Ms i+i(v)£ab(i l x L -i)X a v , where 



s l ( V ) = X l 2 - 1 + (-l) l X[ 



i v i-i 



(A6a) 
(A6b) 

(A7) 

§(i-a), 

±\/l — 4?7 (the ± sign depends on 
one's convention for which mass is larger). 

To compute the "scalar" multipoles defined in Eq. (IA2() 
we need to contract Eqs. (|A6|) with y 1 ™* . Using the re- 
lationship between the "scalar" and STF spherical har- 
monics, 



and we define X\ 



and A = 



M 



|(1+A),X 2 



A/ 



ylrn = ylm nL = ylm ft 



(L), 



(A8) 



the Newtonian "scalar" mass multipole equivalent to 
(|A6a|) is easily seen to be 



iP m =A lV Ms l (r 1 yY 



(A9) 



To derive the Newtonian "scalar" current multipole 
moment we use the definition of the magnetic-type 
"pure-spin" vector spherical harmonics [Eqs. (2.18b) and 
(2.23b) of jH]: 



Y, 



Bin 
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l + l 
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SbcdXc^dYl 
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(A10) 
(All) 



Combining this equation with Eqs. (|A2b|) and (IA6b[) 
yields the Newtonian scalar current multipole for general 
orbits: 



-j-7]Msi +1 (T]y{x x v) ■ VU ( * 



(A12) 



In Eqs. (IA9p and (| A12I) . the moments are given as func- 
tions of time by solving the equations of motion to deter- 
mine the spherical coordinates of the relative separation 
vector x: r(t), 9(t), and (f> = f{t). If we restrict ourselves 
to orbits in the x-y plane, we can further simplify the 
multipole moments by using 



.xx v = v <f>e z , and 
sin 9 dYim 



(A13) 
(A14) 



The resulting Newtonian-order "scalar" multipole mo- 
ments for general orbits restricted to the x-y plane are 



I? m = Air)Msi{rj)r l Y* m (|, <f>) , and 



J iL = -y?7Msi +1 (r?)r 



00 



(A15a) 
(A15b) 



Appendix B: NONLINEAR MEMORY INTEGRAL 
FOR ELLIPTICAL ORBITS IN TERMS OF 
HYPERGEOMETRIC FUNCTIONS 



In this appendix, we show how to derive explicit ex- 
pressions for the integrals of Eqs. (|2.33[) . which we rewrite 
as: 
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where b = f^ff, the Ki m are integration constants, and we refer to the constants in square brackets as Ai m below 
and in the main text. We now note that all of the indefinite integrals in square brackets can be expressed in terms of 



22 



combinations of the following integral [97J: 



„a+l 



— Tudx= — — 2 Fi(b,a+l;a + 2;-cx), 

(1 + cx)° a + 1 

where 2-F1 is the hypergeometric function [98]. Any of the integrals in Eqs. (|Blj) can then be computed from: 
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FIG. 7. (color online). Effect of wavelength averaging on the 
nonlinear memory integrand. The solid (black) curve shows 
the time evolution (parameterized in terms of decreasing e±) 
of the integrand h^ cm ^ 1 ' 1 computed without wavelength av- 
eraging via Eq. (|2.24al) as described in the text. The dashed 
(red) curve shows h^ cm ^ 1 ' 1 computed with wavelength aver- 
aging via Eq. (|2.27a[l . The averaging procedure removes the 
short time scale oscillations from the integrand. The insets 
zoom in on the low and high eccentricity regions. 



where the constants ai r , 



Clr, 



and dim are easily read off 



of Eqs. (|B1[) . and j3 = The integration constants 

Ki m are then determined by the requirement that the 
nonlinear memory vanish at early times when the eccen- 
tricity et 
then given by 



e_. The final result for the /i|™ em ' ) modes is 
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the gravitational- wave energy flux 



dE s 
dtdtt 



This is consis- 



1) _ 



A lm C (e )e* im ™ [F lm (e t ) - F /m (e_)] . (B4) 



tent with the standard derivation in which averaging is 
necessary to obtain a well-defined GW stress-energy ten- 
sor [l6|, OSj. However, in the derivations of the nonlinear 
memory in [bfl . [l5| , this wavelength averaging does not 
explicitly appear. The purpose of this appendix is to in- 
vestigate (in the context of eccentric binaries) how the 
nonlinear memory calculation depends on whether the 
wavelength averaging is performed or not. The short an- 
swer to this question is that the nonlinear memory does 
not depend on this averaging, aside from very small am- 
plitude oscillations at the orbital period that are superim- 
posed on the memory when averaging is not performed. 
The reason why the memory is relatively insensitive to 
the averaging procedure is simple: performing the time 
integral that explicitly appears in the nonlinear mem- 
ory calculation effectively "averages" over the integrand 
[cf. Equation (|2.5[) ]. So by performing also the wavelength 
averaging () of the integrand, one is effectively "averag- 
ing" twice. However, note that the wavelength averaging 
significantly simplifies the integrand, allowing for an an- 
alytic calculationFI 

To investigate this issue in more detail, we can explic- 
itly compute the hi^ onl ^ nonlinear memory mode with 
and without wavelength averaging. In both cases, we 
first solve for the evolution of e t (t), p(t), and v(t) (the 
true anomaly) by numerically integrating Eqs. (|2.29[) and 
(|2. 14c[) . We assume an equal- mass binary and initial con- 
ditions e t (0) = 0.7, p(0) = 30M, and v(0) = 0. We then 
substitute the result into the integrand of the time inte- 
gral for the memory: Eq. (|2.24a[) in the nonaveraged case 
(ignoring the ()), and Eq. (|2.27a[) in the averaged case. 
The resulting integrands are plotted in Fig. [71 There we 
see that if we do not perform any wavelength averaging, 
the integrand retains oscillations at the orbital period; 



Appendix C: THE ROLE OF AVERAGING THE 
GRAVITATIONAL- WAVE STRESS-ENERGY 
TENSOR IN NONLINEAR MEMORY 
CALCULATIONS 

In the definition of the nonlinear memory in Eq. (|2.5I) , 
an explicit averaging over several wavelengths appears in 



Note also that in the circular, nonspinning case, this wavelength 
averaging issue does not arise. The integrand of those modes that 
contribute to the nonlinear memory are constant on an orbital 
time scale and are unaffected by averaging. This can be seen 
explicitly by comparing Eqs. I|2.24|l and l|2.27|l in the et = 
case. We are currently investigating this averaging issue for the 
case of quasicircular, spinning binaries [85ll , 
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FIG. 8. (color online). Effect of wavelength averaging on the nonlinear memory mode. The left-hand plot shows the h.^ " 1 
nonlinear memory mode (the time integral of Fig. [7]). The solid (black) curve is without wavelength averaging; the dashed (red) 
curve is with averaging. The curves lie nearly on top of each other, aside from small amplitude oscillations that remain when 
one integrates the nonaveraged integrand. The inset zooms in on the low-eccentricity region. The integration is terminated at 
the last-stable-orbit. The right-hand plots show the absolute and fractional errors between the two curves. 



these oscillations are smoothed over by the wavelength 
averaging. 

We then numerically integrate both the averaged and 
nonaveraged integrands, starting from the condition 
hf o cm (0) = 0. The result is plotted in the left-hand side of 
Fig. [5) There we see that performing the wavelength av- 
eraging has had very little effect on the resulting memory. 



The two curves lie nearly on top of each other. As stated 
above, this agreement is simply due to the fact that the 
time integration of the nonaveraged /i^™ 6 essentially 
acts as an averaging procedure. The only difference is 
a very small remnant oscillation about the wavelength- 
averaged curve. At the last-stable-orbit and for a large 
range of eccentricity, the two curves agree to ~ 1%. 
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